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1.  INTRODUCTION 


The  penetration  performance  of  high  density,  low  length-to-diameter  ratio  (L/D)  penetrators 
impacting  steel  targets  has  been  a  topic  of  considerable  interest  to  the  armor/antiarmor 
community  due  to  the  speculated  performance  of  a  segmented  rod  penetrator  munition.  This 
has  been  spurred  by  the  observation  that  relative  penetrator  performance  (i.e.,  penetration  per 
unit  length  [P/L])  increases  as  L/D  decreases,  provided  that  impact  velocity  is  relatively  high. 
Most  of  the  research  to  date  which  substantiates  this  claim  has  focused  on  penetrators 
shaped  as  either  spheres  or  right-circular  cylinders  with  an  L/D  of  one  or  slightly  greater 
(Charters,  Menna,  and  Piekutowski  1989;  Charters  1987;  Charters  and  Orphal  1989;  Cline, 
Gogolewski,  and  Reaugh  1987;  Cuadros  1989;  Frank  and  Zook  1990;  Hauver  and  Melani 
1990;  Hohler  and  Stilp  1990;  Hohler  and  Stilp  1987;  Holland  et  al.  1989;  Hunkier  1990;  Kivity, 
Yitzhak,  and  Hirsch  1989;  Naz  and  Lehr  1989;  Raatschen  et  al.  1987;  Scheffler  and  Zukas 
1989;  Scheffler  1989;  Sorensen  et  al.  1991;  Tate  1990;  and  Zukas  1990a).  The  literature 
consists  of  a  mix  of  experiments  and  Computational  simulations  in  which  the  trend  of 
increasing  P/L  with  decreasing  L/D  is  conclusive  for  high  density  metallic  penetrators 
impacting  steel  or  aluminum  targets.  Denardo  (1 989)  finds  this  trend  to  be  true  also  for  low 
L/D  polyethylene  cylinders  impacting  aluminum  targets.  Another  trend  evident  from  the 
literature  is  the  dependence  of  segment  performance  upon  impact  velocity.  Hauver  and 
Melani  (1990);  Kivity,  Yitzhak,  and  Hirsch  (1989);  Orphal,  Anderson,  and  Franzen  (1990); 
Raatschen  et  al.  (1987);  Scheffler  and  Zukas  (1989);  and  Zukas  (1990a)  show  through 
experiments  and  computations  that  low  L/D  segment  performance  increases  dramatically  with 
velocity  and  that  a  train  of  ideally  spaced  segments  does  not  out-perform  its  monolithic  rod 
counterpart  until  the  impact  velocity  approaches  or  exceeds  2  km/s  for  tungsten-RHA  impacts. 
Tate  (1 990)  used  the  modified  Bernoulli  equation  approach  to  develop  an  engineering  model 
for  a  segmented  penetrator  in  which  the  effects  of  segment  spacings  were  also  considered. 

The  performance  of  high  density  metallic  penetrators  with  an  L/D  of  less  than  one  can  be 
estimated  by  extrapolating  the  performance  of  higher  L/D  penetrators.  Unfortunately,  the 
inherent  hazards  of  such  extrapolations  frequently  lead  to  misleading  or  incorrect  results.  A 
limited  number  of  investigations  have  considered  penetrator  shapes  in  which  the  L/D  was  less 
than  one;  however,  these  were  performed  at  very  small  scale.  Herbette  (1 989)  noted  a 
dramatic  increase  in  P/L  for  steel  disks  with  an  L/D  1/30  impacting  aiuminum  targets  at  2  km/s 
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when  compared  to  penetrators  with  considerably  greater  L/D.  Orphal,  Anderson,  and  Franzen 
(1989)  and  Orphal  and  Franzen  (1990)  also  report  a  significant  increase  in  P/L  as  penetrator 
L/D  w P3  reduced  from  1  to  1/8  for  tungsten,  tungsten  alloy,  and  tantalum  alloy  penetrators 
impacting  steel  at  a  velocity  ranging  from  1.5  to  7.5  km/s.  Cnristman  and  Gehring  (1966)  fired 
aluminum  and  steel  penetrators  with  an  L/D  ranging  from  25  to  1/6  into  metallic  and  non- 
metallic  targets;  however,  their  data  was  presented  in  such  a  manner  that  penetrator 
performance  with  relative  L/D  could  not  be  interpreted.  The  Christman  and  Gehring  data  were 
also  subject  to  large,  uncontrolled  penetrator  pitch  and  yaw  at  impact.  Denardo  (1989) 
reported  similar  trends  as  the  L/D  of  polyethylene  penetrators  were  reduced  from  1  to  1/6. 
Despite  these  studies,  a  full  characterization  and  understanding  of  low  L/D  penetrator  behavior 
has  yet  to  be  conducted. 

In  this  study,  a  combined  experimental  and  computational  investigation  was  conducted  to 
determine  the  terminal  ballistic  performance  and  transient  behavior  of  low  L/D  tungsten  alloy 
penetrators  impacting  normal  incidence  rolled  homogeneous  armor  (RHA)  targets  at  2  km/s. 
Each  penetrator  had  a  nominal  diameter  of  25.4  mm,  and  an  L/D  of  either  1, 1/2, 1/4, 1/8, 

1/16,  or  1/32.  Experiments  were  performed  using  a  50-mm  high-pressure  powder  gun  and  a 
push-launch  sabot.  Data  collected  consisted  of  penetrator  attitude  and  velocity  at  impact, 
depth  of  penetration,  crater  diameter,  and  limited  information  regarding  residual  penetrator 
material.  Numerical  simulations  of  each  of  the  penetrators  impacting  RHA  were  performed 
with  the  Lagrangian  code  SMASH  (Kimsey  and  Zukas  1986)  and  the  Eulerian  code  CTH 
(McGlaun  and  Thompson  1990).  Both  codes  exhibited  good  agreement  with  experimental 
data. 


2.  EXPERIMENTS 

A  total  of  13  shots  were  made  in  the  BRL  high  pressure  gun  facility.  Details  of  the  facility 
are  given  by  Bauer  and  Nagy  (1979)  and  Silsby  et  al.  (1983).  The  gun  consists  of  a  nominal 
50-mm-diameter  by  6-m-travel  smooth  bore  powder  gun  with  a  large  capacity,  high  pressure 
powder  chamber.  The  gun  empties  into  an  impact  chamber  at  atmospheric  pressure.  The 
target  was  located  approximately  4.8  m  from  the  gun  muzzle. 
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One  piezoelectric  and  two  copper  crush  pressure  transducers  were  located  in  the  powder 
chamber  of  the  gun  to  aid  with  the  internal  ballistic  diagnostics  of  the  gun.  Two  150-kV 
orthogonal  flash  x-ray  shadowgraph  systems,  located  approximately  450  mm  and  760  mm 
uprange  of  the  target,  were  used  to  provide  penetrator  attitude  and  velocity  data  prior  to 
impact.  The  x-ray  systems  required  a  means  for  detecting  and  enunciating  the  arrival  of  the 
penetrator  at  the  x-ray  stations  (i.e.,  a  properly  timed  trigger  to  activate  each  system).  The 
trigger  system  employed  was  a  "break-screen"  placed  approximately  300  mm  uprange  of  the 
first  x-ray  station.  The  break-screen  consisted  of  a  piece  of  paper  with  a  conductive  etching 
on  it,  mounted  to  a  thin  piece  of  cardboard  located  on  the  shot-line  of  the  gun.  A  voltage  is 
applied  across  the  etching  and  is  monitored  for  interruptions.  When  the  penetrator  strikes  the 
break-screen,  it  perforates  the  cardboard  and  interrupts  the  electrical  circuit  across  the  screen. 
This  interruption  activates  the  delay  generators  used  to  pulse  the  flash  x-ray  systems. 

The  penetrators  were  made  from  a  91%  tungsten  alloy  (Teledyne  Firth  Sterling  Alloy  X-27) 
which  was  swaged  to  25%  reduction  in  area.  The  properties,  static  strengths,  and 
composition  of  the  alloy  are  listed  in  Table  1 . 

The  penetrator  geometries  studied  were  right-circular  cylinders  with  an  L/D  ratio  of  either 
1/2, 1/4, 1/8, 1/16,  or  1/32.  The  diameter  of  each  penetrator  was  nominally  25.4  mm,  except 
for  the  L/D  1/2  penetrator  which  was  29.1  mm.  An  L/D=1  penetrator  with  a  diameter  close  to 
25  mm  was  not  fired  because  it  was  too  massive  for  the  sabot  system  used  for  this  program; 
hence,  its  performance  was  determined  from  calculations.  The  nominal  mass  and  dimensions 
for  each  of  the  different  L/D  penetrators  are  listed  in  Table  2. 

Right-circular,  low  L/D  penetrators  are  highly  unstable  during  flight  and  prone  to  tumbling. 
Conventional  push-  or  traction-launch  sabots  release  the  penetrator  soon  after  muzzle  exit, 
which  provides  ample  time  for  penetrator  tumbling  to  occur,  resulting  in  high  yaw  at  impact.  In 
order  to  achieve  low  yaw  at  impact,  a  different  sabot  technique  was  employed.  A  one-piece 
frangible  sabot  made  from  6/6  Nylon,  Figure  1,  was  used.  The  sabot  had  an  overall  L/D  of 
approximately  one  and  served  as  the  flight  body  for  the  penetrator  after  muzzle  exit.  The 
penetrator  was  placed  in  a  cavity  machined  into  the  front  of  the  sabot.  The  relatively  high- 
mass  penetrator  acted  as  a  nose  weight  which  yielded  an  aerodynamically  stable  sabot- 
penetrator  assembly  .  The  penetrator  was  released  from  the  sabot  on  impact  with  the  break- 
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Table  1.  Tungsten  Alloy  Static  Properties 


Density: 

17,270  kg/m3 

Hardness: 

39.1  HRC 

j  Tensile  Properties 

0.2%  Offset  Yield  Strength: 

1.09  GPa 

!  Ultimate  Strength: 

1.11  GPa 

Total  Elongation  to  Failure: 

9.7% 

1  Compression  Properties  1 

0.2%  Offset  Yield  Strength: 

0.95  GPa 

1 .0%  Offset  Yield  Strength: 

1.16  GPa 

Composition  (by  Weight) 

Tungsten: 

91.0% 

Nickel: 

6.3% 

Iron: 

2.7% 

Table  2.  Penetrator  Geometries 


L/D 

Mass 

(9) 

Thickness 

(mm) 

Diameter 

(mm) 

1/2 

150.4 

12.95 

29.13 

1/4 

55.6 

6.35 

25.40 

1/8 

27.6 

3.17 

25.40 

1/16 

13.8 

1.60 

25.40 

1/32 

6.8 

0.79 

25.40 

screen.  The  impact  forces  scattered  the  Nylon  sabot,  causing  the  debris  to  spread  radially, 
falling  behind  and  leaving  the  penetrator  undisturbed.  The  short  flight  distance  which 
remained  was  not  sufficient  for  penetrator  tumbling  to  occur.  This  technique  provided  a 
reliable  means  for  delivering  disk-shaped  penetrators  to  a  target  with  proper  impact  orientation 
and  negligible  contribution  by  the  sabot  iu  penetration. 

The  targets  were  monolithic  blocks  of  RHA  with  a  nominal  height  and  width  of  200  mm 
and  a  thickness  of  150  mm.  The  measured  hardness  of  the  targets  ranged  from  241  to 
277  BHN.  The  targets  were  placed  in  the  impact  chamber  such  that  the  obliquity  (the  angle 
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SIDE  VIEW 


FRONT  VIEW 


Figure  1.  One-Piece  Frangible  Sabot. 

between  the  target  face  normal  and  the  shot  line)  was  less  than  one  degree  and  the 
penetrator  would  impact  the  center  of  the  target  face. 

The  ballistic  test  procedure  used  was  typical  of  conventional  ballistic  range  impact  testing 
in  that  the  penetrator  was  accelerated  to  the  desired  velocity  in  the  gun,  permitted  to  pass 
through  the  instrumentation  field,  and  its  flight  then  terminated  by  impact  with  the  target. 
Projectile  velocity,  yaw,  and  yaw  rate  measurements  were  obtained  from  the  x-ray 
shadowgraphs  and  associated  timing  information.  The  projectile  velocity  was  assumed  not  to 
decrease  significantly  during  the  short  flight  distance  between  the  x-ray  stations  and  the 
target;  hence,  the  velocity  obtained  from  the  x-ray  systems  was  used  as  the  impact  velocity. 
Penetrator  yaw  at  impact  was  determined  by  assuming  a  constant  yaw  rate  between  the  x-ray 
locations  and  the  target.  A  summary  of  the  experimental  results,  including  velocity  and  yaw, 
is  presented  in  Table  3. 

Post-test  impact  crater  measurements  consisted  of  crater  depth  (P)  and  crater  diameter 
(H).  Impact  crater  depths  were  measured  with  a  depth  micrometer  probe.  The  measured 
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Table  3.  Performance  of  Low  L/D  Penetrators  Impacting  Semi-Infinite  RHA 


Shot 

No. 

L/D 

Velocity 

(m/s) 

Yaw 

(deg) 

P 

(mm) 

P/E 

(mm/MJ) 

P/L 

H 

(mm) 

194 

1/2 

2,042 

4.3 

33.0 

105.2 

66.0 

206 

1/4 

2,100 

4.0 

19.3 

157.4 

49.8 

207 

1/4 

2,044 

1.2 

19.6 

168.7 

3.1 

51.8 

208 

1/4 

2,015 

3.8 

18.1 

160.3 

2.8 

55.9 

216 

1/8 

1,931 

3.0 

13.1 

254.6 

4.1 

37.3 

204 

1/8 

2,040 

16.4 

285.6 

5.2 

39.4 

221 

1/8 

2,036 

18.5 

323.4 

5.8 

41.9 

222 

1/8 

1,996 

13.1 

238.3 

4.1 

41.4 

217 

1/16 

2,004 

10.4 

6.9 

249.0 

29.5 

219 

1/32 

1,958 

12.6 

2.3 

176.5 

27.9 

depths  were  made  relative  to  the  undisturbed  surface  of  the  impacted  target  face.  The 
diameter  of  each  crater  was  also  measured;  however,  this  measurement  was  subjective  due 
to  the  sloping  walls  of  the  relatively  shallow  craters.  The  diameter  was  typically  measured 
slightly  below  the  target  surface.  The  post-test  crater  measurements  are  incluaed  »n  Table  3. 
Penetration  depths  normalized  by  penetrator  length  and  by  penetrator  kinetic  energy  at  impact 
are  also  included  in  Table  3.  A  representative  target  for  each  of  the  different  L/D  penetrators 
(except  L/D  1/2)  was  sectioned  and  photographed  and  are  shown  in  Figure  2.  Residual 
penetrator  material  was  recovered  for  one  shot  with  an  L/D  1/4  penetrator  and  two  shots  with 
an  L/D  1/8  penetrator.  These  were  weighed,  measured,  and  photographed.  The  recorded 
data  are  listed  in  Table  4,  and  a  photograph  of  a  representative  recovered  penetrator  for  L/D 
1/4  and  1/8  is  shown  in  Figure  3.  The  diameter  listed  in  Table  4  is  an  estimated  diameter 
based  upon  the  general  shape  of  the  recovered  penetrator. 


3.  COMPUTATIONS 


Lagrangian  calculations  of  low  L/D  projectiles  impacting  semi-infinite  RHA  were  performed 
with  SMASH,  a  local  variant  of  the  EPIC-2  code  (Johnson  1978)  modified  to  include  eroding 
slide-line  logic  (Kimsey  and  Zukas  1986).  Dynamic  material  properties  for  the  tungsten  alloy 
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L/D  1/8 
SHOT  222 


(b)  Shot  222;  L/D  1/8 

Figure  2.  Sectioned  Target  Crater  Profile 


L/D  1/16 
SHOT  217 


(c)  Shot  217;  L/D  1/16 


(d)  Shot  219;  L/D  1/32 

Figure  2.  Sectioned  Target  Crater  Profiles  (Continued). 
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Table  4.  Residue  Penetrator  Measurements 


Thickness 

(mm) 


Estimated 

Diameter 

(mm) 


Mass 

(g) 


Fraction  of 
Mass  Recovered 


/4 

3 

30 

■■ 

35.4 

/8 

3 

25 

46.4 

18 

3 

25 

17.9 

64.9 

L/D  1/8 


L/D  1/4 


(a)  Shot  204  (b)  Shot  206 

Figure  3.  Residual  Penetrators. 

projectiles  were  taken  from  split-Hopkinson  bar  data  reported  by  Nicholas  (1980)  for 
90W-7Ni-3Fe,  25%  swaged.  The  yield  stress  for  the  tungsten  was  taken  to  be  14  kbar  and 
the  ultimate  stress  to  be  1 9.3  kbar.  The  RHA  target  material  was  modeled  as  an  elastic 
perfectly  plastic  material  with  a  flow  stress  of  7.0  kbar  (Benck  and  Robitaille  1977). 

A  power  series  form  of  the  Mie-Gruneisun  equation-of-state  was  used  to  define  pressure- 
volume  states,  and  the  appropriate  parameters  for  tungsten  alloy  and  steel  were  taken  from 
Kohn  (1969). 
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Two  levels  of  failure  may  be  modeled  with  the  version  of  EPIC  used  for  these  studies, 
based  on  user-specified  levels  of  effective  plastic  strain.  When  the  first  level  is  met,  tensile 
and  shear  stresses  are  not  allowed  to  develop  in  the  affected  elements.  The  net  result  is  that 
the  element  behaves  much  like  a  liquid  in  that  it  can  support  only  hydrodynamic  compression. 
When  the  second  level  of  effective  plastic  strain  is  reached,  ail  stresses  and  pressures  are  set 
to  zero.  Element  quantities  are  no  longer  included  in  the  simulation,  as  though  total  failure 
had  occuned.  When  this  option  is  invoked,  mass  and  momentum  continue  to  be  conserved 
since  the  masses  and  velocities  of  the  failed  elements  are  associated  with  their  nodes,  and 
these  are  tracked  throughout  the  computation.  Energy,  however,  is  conserved  only 
approximately  since  element  internal  energies  are  no  longer  computed.  This  is  a  minor 
concern,  however,  since  this  latter  value  of  strain  is  set  at  high  (150-250%)  levels  (for  lack  of 
better  failure  models)  so  that  elements  "removed"  from  the  calculation  would,  in  reality,  appear 
as  ejecta  in  ballistic  experiments. 

The  grid  for  the  calculations  was  made  up  of  four  constant  strain  triangles  per 
quadrilateral.  The  total  number  of  nodes  for  the  calculation  approached  15,000  while  the  total 
number  of  elements  was  close  to  30,000.  As  a  rule,  between  10  and  25  layers  of 
quadrilateral  elements  were  used  through  the  thickness  of  the  projectile,  the  number  of 
elements  in  the  radial  direction  being  adjusted  to  maintain  a  1:1  aspect  ratio.  The  target  block 
had  a  radius  of  60  mm  and  a  depth  of  100  mm,  effectively  simulating  a  semi-infinite  block. 
Some  calculations  were  redone  with  a  target  radius  of  100  mm,  but  these  showed  no 
appreciable  difference  from  the  calculations  with  a  radius  of  60  mm.  SMASH  input  data  for 
the  cases  modeled  are  presented  in  Appendices  A-F,  and  the  computational  results  are 
shown  in  Table  5. 

A  duplicate  set  of  calculations  were  performed  with  the  Eulerian  code  CTH  (McGlaun  and 
Thompson  1990)  for  comparison.  CTH  is  a  family  of  computer  programs  for  modeling  solid 
dynamics  problems  involving  shock  wave  propagation,  multiple  materials,  and  large 
deformations  in  one,  two,  and  three  dimensions.  The  code  solves  the  finite-difference  analogs 
of  the  conservation  equations  for  mass,  momentum,  and  energy.  A  high  resolution  interface 
tracker  is  used  to  minimize  material  dispersion,  and  a  second-order  convection  scheme  is 
used  to  flux  all  quantities  between  cells.  The  material  model  is  elastic  perfectly  plastic  with 
linear  thermal  softening.  Material  failure  is  based  on  maximum  pressure  or  principal  stress. 
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Table  5.  Computational  Results 


L/D 

L 

(mm) 

P 

SMASH 

(mm) 

P/L 

SMASH 

P 

CTH 

(mm) 

P/L 

CTH 

1 

25.40 

50.8 

2.0 

54.0 

2.1 

1/2 

12.70 

27.9 

2.2 

32.0 

2.5 

1/4 

6.35 

21.6 

3.4 

20.0 

3.2 

1/8 

3.18 

16.2 

5.1 

13.0 

4.1 

1/16 

1.59 

7.5 

4.8 

6.5 

4.1 

1/32 

0.79 

3.0 

3.8 

2.4 

3.0 

Both  analytical  and  tabular  equations-of-state  are  available  for  modeling  the  hydrodynamic 
behavior  of  solid  materials.  In  addition,  a  programmed  bum  model  is  available  of  modeling 
high  explosive  detonation. 


The  CTH  simulations  reported  herein  used  a  linear  Hugoniot  shock-particle  velocity 
equation-of-state  to  model  the  hydrodynamic  behavior  of  the  materials.  An  elastic  perfectly 
plastic  model  is  the  only  plasticity  model  currently  implemented  in  CTH.  Therefore,  the 
dynamic  yield  strengths  for  both  the  penetrator  (19.3  kbar)  and  target  (7.0  kbar)  materials 
used  for  the  Lagrangian  calculations  were  used  for  the  Eulerian  calculations. 


The  CTH  simulations  used  a  cylindrical  coordinate  system  with  a  semi-infinite  RHA  target. 
The  multiple  material  temperatures  and  pressures  thermodynamic  model  was  used  to 
calculate  separate  temperatures  and  pressures  for  each  material  in  multimaterial  cells.  A  high 
resolution  interface  algorithm  was  used  to  define  material  interfaces  and  minimize  material 
dispersion  in  multimaterial  cells. 

The  zoning  for  the  CTH  simulations  employed  square  cells  in  the  penetrator-target 
interaction  region.  Typically,  5-10  cells  were  used  through  the  thickness  of  the  penetrator 
(e.g.,  1.06-mm  by  1.06-mm  cells  were  used  for  the  L/D=1  penetrator  model  and  0.08-mm  by 
0.08-mm  cells  were  used  for  the  L/D=1/32  penetrator  model).  The  total  number  of  cells  in  a 
model  ranged  from  30,000  to  148,125.  Typically,  the  radial  and  axial  boundaries  of  the  target 
were  placed  approximately  5  penetrator  diameters  from  the  impact  point.  Transmitting 
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boundary  conditions  were  used  at  the  radial  and  axial  edges  of  the  mesh.  The  CTH  input 
data  for  all  L/O  cases  modeled  are  provided  in  Appendices  G-L,  and  the  computational  results 
are  summarized  in  Table  5. 

4.  DiSCUSSION 

4.1  Wave  Mechanics:  Correlation  With  Experimental  end  Numerical  Results.  The 
experimental  and  computational  results  indicate  fiat  the  behavior  of  low  L/D  segments  can  be 
loosely  grouped  into  three  categories: 

(1)  rod-like  response  (L/D  >  1) 

(2)  transitional  -  aspects  of  both  rod  and  plate  behavior  (1  >  L/D  >  1/8) 

(3)  plate-like  response  (L/D  <  1/8) 

These  can  best  be  illustrated  using  elementary  concepts  from  wave  mechanics.  For  a 
background  in  this  area,  see  Johnson  (1972),  Wasley  (1973),  Rinehart  (1975),  Zukas  et  al. 
(1982),  and  Zukas  (1990b). 

Consider  a  cylindrical  projectile  with  L/D— 1 .  Upon  impact,  a  compression  wave  is 
propagated  into  both  target  and  projectile  with  an  initial  intensity  of  pcv  where  p  represents 
density,  c  the  sound  speed  in  the  material,  and  v  the  particle  velocity.  For  like-material 
impacts,  v  =  Vs/2,  where  Vs  is  the  striking  velocity.  For  dissimilar  material  impacts,  v  can  be 
determined  from 

v  =  Vs  /  [1  +  ( pt  c, )  /  (  pp  cp )] , 

where  subscripts  t  and  p  refer  to  target  and  projectile,  respectively.  For  a  tungsten  alloy 
projectile  (p  =  17.3  g/cm3;  c  =  0.4  cm/ps)  striking  an  RHA  plate  (p  =  7.8  g/cm3; 
c  =  0.55  cm/ps)  at  a  velocity  of  2  km/s,  v  =  0.12  cm/ps,  and  the  pressure  in  the 
penetrator  =  pcv  =  830  kbar.  Almost  immediately  upon  impact,  relief  waves  propagate  in  from 
the  lateral  free  surfaces.  After  a  time  t  =  r/c  »  3.2  ps,  relief  waves  from  the  radial  boundary 
will  have  reached  the  center  of  the  projectile  and  diminished  the  magnitude  of  the  initial 
compressive  pulse.  The  original  compressive  pulse  will  reach  the  rear  surface  of  the  projectile 
at  a  time  t  =  L/c  »  6.4  ps.  By  this  time,  a  fully  two-dimensional  stress  state  will  have  been 
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developed  in  the  projectile  with  a  pressure  amplitude  considerably  less  than  the  peak 
pressure  first  experienced  on  impact.  This  is  shown  in  Figures  4  and  5,  which  are  plots  of 
pressure  as  a  function  of  time  at  two  locations  in  the  projectile— the  projectile-target  interface 
and  the  rear  surface  of  the  rod. 

Shortly  aner  t  =  L/c,  the  pulse  is  reflected  from  the  rear  of  the  projectile  and  returns  to  the 
projectile-target  interface  at  t  =  2L/c.  Thence,  due  to  the  impedance  mismatch  between  the 
projectile  and  target,  a  portion  of  the  pulse  is  transmitted  into  the  target,  the  remainder  being 
reflected  back  into  the  projectile.  Wave  reverberations  in  the  radial  direction  in  the  projectile 
continue  as  well  so  that,  until  consumed  or  arrested,  the  L/D=1  projectile  will  experience  a 
two-dimensional  stress  state  which,  after  the  first  wave  reverberations,  is  too  complex  to 
describe  analytically  but  can  readily  be  tracked  numerically  (e.g.,  Zukas  and  Sf/gletes  1989). 

Note  that  the  computer  calculations  account  for  two-dimensional  geometric  effects, 
pressure  pulses  with  finite  rise  times,  and  changes  in  sound  speed  with  compression.  Thus, 
there  will  be  small  differences  in  arrival  times  between  those  from  elementary  wave  theory 
and  numerical  computations.  The  basic  phenomena,  however,  remain  the  same. 

As  the  projectile  length  and  L/D  ratio  decreases  (D  being  fixed  in  this  case),  there  begins  a 
transition  from  rod-like  behavior  (two-dimensional  stress  state)  which  is  characteristic  of  long 
rod  penetrators  to  piate-like  behavior  (uniaxial  strain)  for  a  significant  portion  of  the  plate- 
projectile  response.  This  can  be  appreciated  from  Table  6,  which  compares  the  number  of 
wave  transits  <n  the  radial  direction  (r  =  1 .27  mm  for  all  cases)  to  the  number  of  wave  transits 
in  the  axial  direction  with  decreasing  l_/D.  Thus,  for  L/D=1,  in  the  time  the  radial  release  wave 
has  traversed  a  distance  r,  the  axial  wave  has  only  completed  1/4  of  its  passage  through  its 
path  2L.  By  contrast,  for  L/D=  1/32,  there  will  be  eight  complete  transits  in  the  axial  direction 
before  the  arrival  at  the  projectile  center  of  the  radial  release  wave. 

Figure  6  shows  the  pressure-time  history  at  the  rear  of  the  L/D=1/2  and  1/4  projectiles.  A 
large  compressive  pulse  arrives  at  the  rear  surface  and  is  immediately  reflected  in  tension. 

This  tensile  state  is  maintained  for  a  considerable  period,  being  reinforced  by  tension  from  the 
radial  release  wave,  resulting  in  spallation  of  the  projectiles  (illustrated  in  Figure  7).  Residual 


13 


•jQ  2-2  4  0  6-.'  2  0  n,0  12  0  14  0  ISO  15  0  20  0 
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(b)  L/D  =  1/4 


Figure  6.  Predicted  (SMASH)  Pressure  vs.  Time  History  for  Rear  Surface  of  UD=V2  and  1/4 
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Table  6.  Wave  Transits  in  Cylindrical  Projectiles 


L/D 

No.  of  Radial 
Transits 

No.  of  Axial 
Transits 

1/1 

1 

1/4 

1/2 

1 

1/2 

1/4 

1 

1 

1/8 

1 

2 

1/16 

1 

4 

1/32 

1 

8 

Note:  Distance  traveled  in  the  radial  direction  -  r  -  1.27  mm. 
Distance  traveled  in  the  axial  direction  «  2L. 


projectiles  were  recovered  for  the  L/D=1/4  and  1/8  experiments  and  measured  0.3  cm  thick, 
compared  to  numerically  computed  thicknesses  of  0.25  cm. 

For  L/D  £  1/8,  a  state  of  uniaxial  strain  is  approached  in  the  central  region  of  the  projectile. 
The  decreased  projectile  length  permits  numerous  wave  reverberations  in  the  axial  direction 
before  the  arrival  of  the  radial  wave  so  that  an  effective  uniaxial  strain  state  is  maintained  for 
quite  some  time.  Rapid  unloading  of  the  initial  peak  pressure  occurs  (see  Figure  8),  and  the 
short  duration  of  high  amplitude  pressure  waves  inhibit  spall.  There  is  slight  radial  expansion 
of  the  projectile  but  little  to  no  contraction  in  the  axial  direction.  This  is  evident  from  a  time 
history  plot  of  the  front  and  rear  surfaces  of  L/D  1/16  and  1/32  penetrators  (shown  in 
Figures  9  and  10).  The  vertical  distance  between  the  front  and  rear  surface  curves  on  each 
plot  indicate  penetrator  thickness.  The  penetration  crater,  which  for  the  L/D=1  case  resembles 
that  of  a  stubby  cylindrical  projectile  (Figure  11),  is  flat  for  the  L/D  1/32  case  (Figure  12). 
Considerable  penetrator  flexure  was  also  observed  for  the  L/D  <  1/8  numerical  simulations. 

An  example  of  the  mode  of  flexure  is  shown  in  Figure  12a.  A  section  of  the  penetrator  near 
the  radial  center  temporarily  looses  contact  with  the  target  surface.  This  was  observed  in  both 
codes  but  was  more  pronounced  in  the  Lagrangian  simulations. 

The  experimental  and  computational  results  show  that  the  maximum  penetration  per  unit 
length,  or  length  efficiency  of  the  penetrator  (P/L),  is  achieved  for  penetrators  with  an  L/D  of 
about  1/8.  A  plot  of  P/L  values  obtained  from  experiment  and  computations  as  a  function  of 
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Figure  9.  Front  and  Rear  Surface  Displacement  L/D  1/16  Projectile. 


18 


19 


(a)  SMASH  Results 
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Figure  12.  Computed  Penetration  Crater.  LTD  1/32  Fenetrator. 
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L/D  is  sho.vn  in  Figure  13.  A  moderate  degree  of  scatter  exists  in  the  experimental  data, 
especially  for  the  L/D  1/8  shots.  There  is  only  a  slight  discrepancy  in  the  results  obtained 
from  the  two  computer  codes.  However,  experiments  and  computations  show  the  same 
trends  and  agree  quite  well. 

Another  way  of  measuring  penetrator  "efficiency"  is  to  consider  penetration  normalized  by 
penetrator  kinetic  energy,  P/E,  at  impact.  P/E  is  related  to  P/L  by  the  following  relation: 

P/E  =  P/L  [8  /  ( %  pp  D2  Vs2 )] 

Penetrator  diameter  and  impact  velocity  were  held  constant  in  this  study;  hence,  penetrator 
energy  efficiency  is  linearly  related  to  length  efficiency.  Therefore,  Figure  13  also  indicates 
the  trend  of  P/E  with  velocity.  As  can  be  seen  in  Figure  13  (and  Table  3),  a  penetrator  with 
an  L/D  of  1/32  is  more  energy  efficient  than  a  penetrator  with  an  L/D  of  1/2  or  1,  and  the 
penetrator  with  the  greatest  energy  efficiency  is  one  with  an  L/D  of  approximately  1/8. 

An  important  point  to  note  is  that  the  ratio  P/L  does  NOT  tend  to  infinity  as  L  approaches 
zero.  The  transition  from  rod  mechanics  to  plate  mechanics  should  preclude  this,  at  least  at 
velocities  at  or  below  2  km/c.  The  plot  of  P/L  vs  L/D  (Figure  13)  shows  a  peak  of  non- 
dimensional  penetration  at  L/D=1/8,  then  a  decline.  This  suggests  that  while  penetration 
tends  to  zero  as  L  tends  to  zero,  there  is  a  point,  dependent  upon  projectile  material, 
geometry,  and  striking  velocity,  where  P  approaches  zero  at  a  greater  rate  than  L  approaches 
zero.  A  similar  argument  can  be  made  on  the  basis  of  projectile  kinetic  energy.  As  projectile 
mass  tends  to  zero,  projectile  energy  tends  to  zero.  Thus,  projectile  penetration  must 
inevitably  tend  to  zero.  The  case  of  zero  projectile  mass  will  never  be  achieved  in  practice. 
Nevertheless,  with  reduced  projectile  mass  (decreasing  L),  the  accompanying  penetration  will 
decline  from  a  measurable  quantity  to  mere  scratches  on  the  target  surface. 

4.2  Numerical  Considerations.  As  Figures  4,  5,  and  8  indicate,  the  transient  loading 
which  results  from  the  high-velocity  impact  of  a  low  L/D  tungsten  alloy  projectile  with  an  RHA 
target  is  analogous  to  shock  loading.  There  is  a  sharp  compressive  pressure  spike  followed 
by  rapid  tensile  unloading  (both  events  occurring  in  a  time  frame  of  2L/c)  and  then  oscillations 
about  a  mean  value  near  the  material  strength  which  last  for  the  bulk  of  the  response  time. 
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Figure  13.  Experimental  and  Computational  P/L  and  P7E  vs.  L/D  for  2  km/s  Impact  Velocity. 

Adequate  resolution  of  this  pressure  spike  requires  a  very  fine  spatial  grid  in  the  projectile. 
Most  calculations  with  SMASH  were  run  with  at  least  10  elements  through  the  projectile 
length,  the  number  of  elements  in  the  radial  direction  being  adjusted  to  maintain  a  1 :1  aspect 
ratio.  With  this  spatial  resolution,  SMASH  predicted  a  peak  impact  pressure  of  850  kbar, 
compared  to  830  kbar  predicted  from  elementary  wave  theoiy.  CTH  computed  a  peak  value 
of  735  kbar.  Considering  the  differences  between  the  Lagrange  and  Euler  computational 
techniques  as  well  as  the  different  forms  of  the  Mie-GrOneisen  equation-of-state,  the  results 
from  both  codes  are  equivalent. 

Failure  to  maintain  adequate  spatial  resolution  has  other  deleterious  effects.  In  Euler 
calculations,  diffusion  increases  with  increasing  cell  size,  leading  to  the  presence  of  more  and 
more  mixed  cells.  While  the  multiple  material  temperature  and  pressure  thermodynamic 
model  was  used  in  the  Euler  calculations,  pressure-time  histories  at  user-specified  data- 
collection  points  were  obtained  from  interpolations  between  cell-centered  values  which  only 
approximate  the  pressure  at  material  interfaces.  Lagrange  calculations  employing  eroding 
slideline  logic,  such  as  SMASH,  can  suffer  from  excessive  material  loss  and  pressure 
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oscillations  at  the  interface.  For  calculations  which  use  a  relatively  coarse  grid,  the  pressure¬ 
time  history  is  computed  incorrectly,  and  the  mass  removed  by  the  erosion  algorithm  is 
overestimated.  The  zoning  used  in  projectile  and  target  in  these  numerical  simulations  is 
believed  to  be  the  minimum  required  to  accurately  predict  the  residual  length  or  spall 
thickness  as  well  as  peak  impact  pressure. 

The  key  to  successful  computations  of  the  impact  response  of  low  L/D  projectiles  at  high 
velocity  is 

•  high  spatial  resolution 

•  material  data  for  the  constitutive  law  from  dynamic  (wave  propagation)  experiments. 

The  former  implies  a  minimum  of  10  elements  through  the  projectile  length  and  maintaining  a 
1 :1  aspect  ratio  of  the  elements  in  zones  undergoing  large  deformations  or  experiencing  steep 
pressure  gradients.  To  minimize  spurious  pressure  oscillations  resulting  from  the  contact- 
impact  (slideline)  algorithm,  approximately  equal  nodal  masses  should  be  maintained  on  each 
side  of  a  slideline.  This  requirement  imposes  a  restriction  on  the  maxim  m  element  size  in 
the  target  (assuming  the  projectile  grid  to  be  fixed)  in  the  vicinity  of  steep  gradients  of 
pressure  or  deformation  and,  thus,  sets  the  requirement  for  the  total  number  of  elements  in  a 
calculation.  For  good  agreement  with  experiment,  nearly  30,000  elements  were  required  for 
the  Lagrangian  calculations  in  both  projectile  and  target.  For  most  calculations  fewer  than 
2,000  elements  were  required  to  describe  projectile  geometry.  The  Eulerian  calculations  used 
between  30,000  and  148,125.  The  requirement  for  dynamic  material  property  data  for  the 
penetrator  and  target  materials  was  satisfied  by  using  split-Hopkinson  bar  data,  determined  at 
strain  rates  between  100  and  I.OuO'1  s  (Nicholas  1980)  for  the  deviatoric  behavior,  and  a 
high-pressure  equation-of-state  based  on  Hugoniot  data  for  the  volumetric  response  (Kohn 
1969). 

5.  CONCLUSIONS 

For  the  cases  studied  here  (impact  velocity  of  2  km/s;  tungsten  alloy  projectiles  vs.  semi¬ 
infinite  RHA;  projectile  length-to-diameter  ratios  ranging  from  1  to  1/32),  the  following 
conclusions  apply: 
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(1)  The  mechanics  of  deformation  for  the  projectile  change  from  rod-like  behavior  to  plate¬ 
like  behavior  as  L/D  decreases. 

(2)  Projectiles  with  an  L/D  =  1/2  and  1/4  are  observed  to  spall.  For  L/D  £  1/8,  the 
pressure  is  attenuated  too  rapidly  for  spall  to  occur.  For  projectiles  with  L/D  £  1,  relief 
waves  from  the  lateral  free  surfaces  provide  stress  relief  which  eliminate  projectile 
spalling. 

(3)  Projectiles  with  L/D  £  1/8  show  slight  lateral  expansion  and  negligible  length  reduction. 

(4)  Projectiles  with  L/D  £  1/8  show  considerable  flexure  in  which  the  central  radial  region 
of  the  penetrator  temporarily  looses  contact  with  the  target  surface  during  penetration. 

(5)  Absolute  penetration  decreases  with  decreasing  L/D  for  constant  diameter  penetrators. 

(6)  P/L  increases  with  decreasing  L/D,  BUT  ONLY  UP  TO  A  POINT  (in  this  case,  to  an 
L/D  ratio  of  approximately  1/8).  Beyond  that,  there  is  a  decrease  in  P/L  with 
decreasing  L/D.  This  decrease  can  be  related  to  the  transition  from  rod  mechanics  to 
plate  mechanics  in  the  projectile. 
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31 


Intentionally  left  blank. 


32 


120  0  1  1  15000. 

7.8  17.6  2.785 

1.00E+20  1.00E+20  1.00E+20 

7.78E+1 1  1.28E+12  2.60E+11 

0.0  0.0  0.0 

7.00E+09  14.3E+09  2.9E+09 

7.00E+09  19.0E+09  2.9E+09 

0.01  0.10  0.40 

3.50E+10  2.00E+10  2.00E+10 

1.019E+12  3.021  £+12  0.76762E+12 

4.220E+12  4.699E+12  1.28312E+12 

5.129E+12  3.350E+12  1.25133E+12 


1.69 

1.43 

2.00 

0.10 

0.10 

0.10 

2.00 

2.00 

2.00 

WOO 

1.00 

1.0 

1.00 

1.00 

1.00 

2.2 

2.2 

2.20 

21. 

21. 

21. 

1 

0  0 

0 

1.0 

1.0 

0. 

0.0  0. 

6 

1  14 

27  2.54 

0.0 

1.0  1.0  0.0  0.0  0.0 
6  1001  61  101  0.0  -10.  6.00 

6  1  13  26  1  2 


0.0 

0.0 


1.0 


1.000 


1.0 


1.0 


6  1001  60  100  2001  1 

40  22  -1  2 

2  1 

14  742  28 
741  729  -1 
1001  1022  1 
0.0  -2.00E+05 


1.0E-10  1.0E-07  1.0E-11  0.6  80.E-06 


0.E-06 

1111. 

1111. 

0.0 

1 

5.0E-06 

1111. 

1111. 

0.0 

4 

f 

10.E-06 

1111. 

1111. 

0.0 

1 

15.E-06 

1111. 

1111. 

0.0 

1 

20.E-06 

1111. 

1111. 

0.0 

? 

30.E-06 

1111. 

1111. 

0.0 

1 

40.E-06 

1111. 

1111. 

0.0 

1 

50.E-06 

1111. 

1111. 

0.0 

1 

60.E-06 

1111. 

1111. 

0.0 

1 

80.E-06 

1111. 

1111. 

0.0 

1 

800.E-06 

1111. 

1111. 

0.0 

1 

900.E-06 

1111. 

1111. 

0.0 

1 

1000.E-06 

1111. 

1111. 

0.0 

1 

1100.E-06 

1111. 

1111. 

0.0 

1 

1200.E-06 

1111. 

1111. 

0.0 

1 

1400.E-06 

1111. 

1111. 

0.0 

1 

1600.E-06 

1111. 

1111. 

0.0 

1 

1.0 


34 


APPENDIX  B: 
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L/D  1/4  SMASH  INPUT  DECK 
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0.0 

42 


APPENDIX  D: 

L7D  1/8  SMASH  INPUT  DECK 


43 


Intentionally  left  blank. 


44 


120  0  1  1  15000. 

7.8  17.6  2.785 

1.00E+20  1.00F+20  1.00E+20 

7.78E+1 1  1 .28E+1 2  2.60E+1 1 

0.0  0.0  0.0 
7.00E+C9  14.3E+09  2.9E+09 

7.00E+09  19.0E+09  2.9E+09 

0.01  0.10  0.40 

3.50E+10  2.00E+10  2.00E+10 

1.019E+12  3.021  E+1 2  0.76762E+12 

4.220E+12  4.699E+12  1.28312E+12 

5.129E+12  3.350E+12  1.25133E+12 


1.69 

1.43 

2.00 

0.10 

0.01 

0.10 

2.00 

2.00 

2.00 

1.00 

1.00 

1.0 

1.00 

1.00 

1.00 

2.2 

2.2 

2.20 

21. 

21. 

21. 

1 

0  0 

0 

1.0 

1.0 

0. 

0.0  0.00 

6 

1  41 

11  .3175 

0.0  1.270 

1.0  1.0  0.0  0.0  0.G 

6  1001  61  101  0.0  -10.  6.00  0.0 

6  1  40  10  12 


1.0 


1.000 


1.0 


1.0 


6  1001  60  ICO  2001  1 

51  22  -1  2 
2  1 

41  861  82 
860  821  -1 
1001  1022  1 
0.0  -2.00E+05 


1.0E-10  1.0E-07  1.0E-11  0.6  80.E-06 


O.E-06 

1111. 

1111. 

0.0 

1 

5.0E-06 

1111. 

1111. 

0.0 

1 

1 0.E-06 

1111. 

1111. 

0.0 

1 

15.E-06 

1111. 

1111. 

0.0 

1 

20.E-06 

1111. 

1111. 

0.0 

1 

30.E-06 

1111. 

1111. 

0.0 

1 

40.E-06 

1111. 

1111. 

0.0 

1 

50.E-06 

1111. 

1111. 

0.0 

1 

60.E-06 

1111. 

1111. 

0.0 

1 

80.E-06 

1111. 

1111. 

0.0 

1 

800.E-06 

1111. 

1111. 

0.0 

1 

900.E-06 

1111. 

1111. 

0.0 

1 

1000.E-06 

1111. 

1111. 

0.0 

1 
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1100.E-06 

1111. 

1111. 

0.0 

1200.E-06 

1111. 

1111. 

0.0 

1400.E-06 

1111. 

1111. 

0.0 

1600.E-06 

1.0 

1111. 

1111. 

0.0 

46 


APPENDIX  E: 

IVD  1/16  SMASH  INPUT  DECK 


47 


Intentionally  left  blank. 


48 


120  0  1  1  72000. 

7.8  17.6  2.785 

1.00E+20  1.00E+20  1.00E+20 

7.78E+1 1  1.28E+12  2.60E+11 

0.0  0.0  0.0 

7.00E+09  14.3E+09  2.9E+09 

7.00E+09  19.0E+09  2.9E+09 

0.01  0.10  0.40 

3.50E+1 0  2.00E+10  2.00E+10 

1.019E+12  3.021  E+1 2  0.76762E+12 

4.220E+12  4.699E+12  1.28312E+12 

5.129E+12  3.350E+12  1.25133E+12 


1.69 

1.43 

2.00 

0.10 

0.10 

0.10 

2.00 

2.00 

2.00 

989. 

999. 

1.0 

1.00 

1.00 

1.00 

2.2 

2.2 

2.20 

21. 

21. 

21. 

10  0  0 

1.0  1.0  0.  0.0  0.00 

6  1  85  11  .15875  0.0  1.270  0.0  1.0  1.0 

1.0  1.0  0.0  0.0  0.0 

6  2001  79  79  0.0  -6,0  6.00  0.0  1.000  1.0 


6  1  84  10  12 


6  2001  78  78  3401  1 

95  22  -1  2 
2  1 

85  1785  170 
1784  1701  -1 
2001  2022  1 
0.0  -2.00E+05 


1.0E-10 

1.0E-07  1.0E-11  0.6 

100.E-06 

0.E-06 

1111. 

1111. 

0.0 

1 

5.0E-06 

1111. 

1111. 

0.0 

1 

10.E-06 

1111. 

1111. 

0.0 

1 

15.E-06 

1111. 

1111. 

0.0 

1 

20.E-06 

1111. 

1111. 

0.0 

1 

25.E-06 

1111. 

1111. 

0.0 

1 

30.E-06 

1111. 

1111. 

0.0 

1 

40.E-06 

1111. 

1111. 

0.0 

1 

100.E-06 

1111. 

1111. 

0.0 

1 

200.E-06 

1111. 

1111. 

0.0 

1 

300.E-06 

1111. 

1111. 

0.0 

1 

500.E-06 

1111. 

1111. 

0.0 

1 

1000.E-06 

1111. 

1111. 

0.0 

1 
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1100.E-06 

1111. 

1111. 

0.0 

1200.E-06 

1111. 

1111. 

0.0 

1400.E-06 

1111. 

1111. 

0.0 

1600.E-06 

1.0 

1111. 

1111. 

0.0 

50 


APPENDIX  F: 

L/D  1/32  SMASH  INPUT  DECK 


51 


Intentionally  left  blank. 


52 


120  0  1  1  72000. 

7.8  17.6  2.785 

1.00E+20  1.00E+20  1.00E+20 

7.78E+1 1  1.28E+12  2.60E+11 

0.0  0.0  0.0 

7.00E+09  14.3E+09  2.9E+09 

7.00E+09  19.0E+09  2.9E+09 

0.01  0.10  0.40 

3.50E+10  2.00E+10  2.00E+10 

1.019E+12  3.021  E+1 2  0.76762E+12 

4.220E+12  4.699E+12  1.28312E+12 

5.129E+12  3.350E+12  1.25133E+12 


1.69 

1.43 

2.00 

0.10 

0.10 

0.10 

2.00 

2.00 

2.00 

999. 

999. 

1.0 

1.00 

1.00 

1.00 

2.2 

2.2 

2.20 

21. 

21. 

21. 

10  0  0 

1.0  1.0  0.  0.0  0.00 

6  1  85  11  .079375  0.0  1.270  0.0  1.0  1.0 

1.0  1.0  0.0  0.0  0  0 

6  2001  61  101  0.0  -10.  10.0  0  0  1.000  1.0 


6  1  84  10  12 


6  2001  60  100  3401  1 

95  22  -1  2 

2  1 

85  1785  170 
1784  1701  -1 
2001  2022  1 
0.0  -2.00E+05 

1.0E-10  1.0E-07  1.0E-11  0.6  100.E-06 


0.E-06 

1111. 

1111. 

0.0 

1 

5.0E-06 

1111. 

1111. 

0.0 

1 

10.E-06 

1111. 

1111. 

0.0 

1 

15.E-06 

1111. 

1111. 

0.0 

1 

20.E-06 

1111. 

1111. 

0.0 

1 

25.E-06 

1111. 

1111. 

0.0 

1 

30.E-06 

1111. 

1111. 

0.0 

1 

40.E-06 

1111. 

1111. 

0.0 

1 

100.E-06 

1111. 

1111. 

0.0 

1 

200.E-06 

1111. 

1111. 

0.0 

1 

300.E-06 

1111. 

1111. 

0.0 

1 

500.E-06 

1111. 

1111. 

0.0 

1 

1000.E-06 

1111. 

1111. 

0.0 

1 
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1100.E-06 

1111. 

1111. 

0.0 

1 

1200.E-06 

1111. 

1111. 

0.0 

1 

1400.E-06 

1111. 

1111. 

0.0 

1 

1600.E-06 

1111. 

1111. 

0.0 

1 

1.0 
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APPENDIX  G: 

L/D  1  CTH  INPUT  DECK 


55 


Intentionally  left  blank. 


56 


*eor*  cgenin 
* 

cth  l/d=1  wa  impacting  rha  at  2km/s 
★ 

control 

ep 

mmp 

*  usessd 
endcontrol 
convct 

convection  1 
interface  high 

endc 

* 

mesh 

block  1  geom  2dc  type  e 
xO  0. 

xl  n  150  dxf  0.1058332  rat  1.0 
endx 

yO  *3.598332 

yl  n  200  dyf=.1 058332  rat  1.0 
endy 

xact  0.0  2.0 

yact  -3.0  1 .0 

endblock 

endmesh 

* 

insertion  of  material 
block  1 

package  ’target' 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yl  0.0  y2  1000.0 

endinsert 

endpackage 

package  ’rod’ 

material  2 

numsub  25 

velocities  yvel  2.Ge5 

insert  box 

xl  0.0  x2 1.27 

yl  -2.54  y2  0.0 

endinsert 

endpackage 

endblock 

endinsert  ion 
* 

edit 
block  1 
expanded 


endblock 

ended  it 
* 

eos  num  2 
matl  -1 
mat2  -2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1.69  0. 
aneosi  3.1416  0.  -1.  -7.36e10  0.  0.  0. 
aneosi  0.  0.  0.  0.  0.  0.  0. 
aneosi  26 1. 

aneosi  -2  tungsten  rhug=-1  thug=-1 

aneos2  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneos3  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneos4  0.  0.  0.  0.  0.  0.  0. 

aneos5  74.  .9 

aneos6  28  .07 

aneos7  26  .03 

endeos 

* 

tracer 

block  1 

add  0.0  -0.15 

add  0.0  -2.35 

endtracer 

* 

endinput 

* 

*eor*  cthin 

cth  l/d=1  wa  impacting  rha  at  2km/s 
* 

control 

mmp 

tstop  10.0e-05 
cpshift  999. 
rdumpf=14400. 

*  pvoid  10.0 

*  usessd 

endcontrol 

* 

edit 

shortt 

time  0.  dtfrequency  1 .0e-06 
ends 
longt 

time  0.  dtfrequency  1 .0e-05 
endl 
plott 

time  0.  dtfrequency  1 .0e-05 
endp 
histt 

time  0.  dtfrequency  1 .0e-07 
time  1.0e-05  dtfrequency  1 .0e-06 


htracerl 

htracer2 

endh 

?nde 

* 

boundary 

bhydno 

block  1 

bxbot  0 

bxtop  1 

bybot 1 

bytop  1 

endb 

endh 

endb 

* 

epdata 

matep  1  yield  7.0e9  poisson  0.33 
matep  2  yield  19.3e9  poisson  0.30 
mix  3 

ende 

* 

tracts 
pressure 
pfmix  -30.0e9 
pfvoid  -30.0e9 
pfracl  -35.0e9 
pfrac2  -20.0e9 
endf 


Intentionally  left  blank. 
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APPENDIX  H: 

L/D  1/2  CTH  INPUT  DECK 


61 


Intentionally  left  blank. 


62 


* 

*eor*  cgenin 
* 

cth  i/d=.50  wa  impacting  rha  at  2knVs 
* 

control 

ep 

mmp 

*  usessd 
endcontrol 
convct 

convection  1 
interface  high 

endc 

* 

mesh 

block  1  geom  2dc  type  e 
xO  0. 

xl  n  300  dxf  0.0529166  rat  1.0 
endx 

yO -3.069166 

yl  n  224  dyf=.0529l66  rat  1.0 
y2  n  20  dyf=.0529166  rat  1.05 
endy 

xact  0.0  2.0 

yact  -1.5  1.0 

endblock 

endmesh 

* 

insertion  of  material 
block  1 

package  'target' 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yl  0.0  y2  1000.0 

endinsert 

endpackage 

package  'rod' 

material  2 

numsub  25 

velocities  yvel  2.0e5 

insert  box 

xl  0.0  x2  1.27 

yl  -1.27  y2  0.0 

endinsert 

endpackage 

endblock 

endinsert  ion 
•* 

edit 
block  1 
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expanded 

endbtock 

endedit 

* 

eos  num  2 
mail  *1 
mat2  -2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1.69  0. 
aneosi  3.1416  0.  -1.  -7.36e10  0.  0.  0. 
aneosi  0.  0.  0. 0.  0.  0.  0. 
aneosi  26  1 . 

aneosi  -2  tungsten  rhug=-1  thug=-1 

anoos2  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneos3  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneos4  0.  0.  0.  0.  0.  0.  0. 

aneosS  74.  .9 

aneos6  28  .07 

aneos7  26  .03 

endeos 

* 

tracer 

block  1 

add  0.0  -0.05 

add  C.O  -1.26 

endtracer 

* 

endinput 

* 

*eor*  cthin 

cth  l/d=.5  wa  impacting  rha  at  2km/s 

<r 

control 

mmp 

tstop  10.0e-05 
cpshift  999. 
rdumpf=14400. 

*  pvoid  10.0 

*  usessd 
endcontrol 

t, 

edit 

shorn 

time  0.  dtfrequency  1 .0e-06 
ends 
longt 

time  0.  dtfrequency  1 .0e-05 
endl 
plott 

time  0.  dtfrequency  1 .0e-C5 
endp 
histt 

time  0.  dtfrequency  1.0e-07 


timel.Oe-05  dtfrequency  1 .0e-06 

htracerl 

htracer2 

endh 

ende 

* 

boundary 

bhydro 

block  1 

bxbot  0 

bxtop  1 

bybot 1 

bytop  1 

endb 

endh 

endb 

* 

epdata 

matep  1  yield  7.0e9  poisson  0.33 

matep  2  yield  19.3e9  poisson  0.30 

mix  3 

ende 

* 

tracts 

pressure 

pfmix  -30.0e9 

pfvoid  -30.0e9 

pfracl  -35.0e9 

pfrac2  -20.0e9 

endf 

* 
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Intentionally  left  blank. 


66 


APPENDIX  I: 

L/D  1/4  CTH  INPUT  DECK 


67 


Intentionally  left  blank. 


68 


*eor*  cgenin 
* 

cth  l/d=.250  wa  impacting  rha  at  2km/s 
* 

control 

ep 

mmp 

*  usessd 
endcontrol 
convct 

convection  1 
interface  high 

endc 

* 

mesh 

block  1  geom  2dc  type  e 
xO  0. 

xl  n  300  dxf  0.0529166  rat  1.0 
endx 

yO  -1.2699984 

yl  n224dyf=.0529166  rat  1.0 
endy 

xact  0.0  2.0 

yact  -1.0  1.0 

endblock 

endmesh 

* 

insertion  of  material 
block  1 

package  'target' 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yl  0.0  y2  1000.0 

endinsert 

endpackage 

package  'rod' 

material  2 

numsub  25 

velocities  yvel  2.0e5 

insert  box 

xl  0.0  x2 1.27 

yl  -.635  y2  0.0 

endinsert 

endpackage 

endblock 

endinsertion 

* 

edit 
block  1 
expanded 


endbtock 

endedit 

* 

eos  num  2 
matl  -1 
mat2-2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1.69  0. 
aneosi  3.1416  0.  -1.  -7.36^10  0.  0.  0. 
aneosi  0.  0. 0.  0.  0.  0.  0. 
aneosi  26  1. 

aneosi  -2  tungsten  rhug=-1  thug=-1 

aneos2  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneos3  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneos4  0.  0.  0.  0.  0.  0.  0. 

anecs5  74.  .9 

aneos6  28  .07 

aneo$7  26  .03 

endeos 

* 

tracer 

block  1 

add  0.0  -0.01 

add  0.0  -0.634 

add  0.0  -0.3175 

endtracer 

* 

endinput 

* 

*eor*  cthin 

cth  l/d*.25  wa  impacting  rha  at  2km/s 
* 

control 

mmp 

tstop  10.0e-05 
cpshift  999. 
rdumpf=14400. 

*  pvoid  10.0 

*  usessd 

endcontrol 

* 

edit 

shortt 

time  0.  dtfrequency  1 ,0e-06 
ends 
longt 

time  0.  dtfrequency  1 ,0e-05 
endl 
plott 

time  0.  dtfrequency  1 .0e-05 
endp 
histt 
time  0. 


dtfrequency  0.05e-06 


time  2.0e-06  dtfrequency  0.1e-06 

time  IO.e-06  dtfrequency  1 .0e-06 

htracerl 

htracer2 

htracer3 

endh 

ende 

* 

boundary 

bhydro 

block  1 

bxbot 0 

bxtop  1 

bybot 1 

bytop  1 

endb 

endh 

endb 

* 

epdata 

matep  1  yield  7.0e9  poisson  0.33 

matep  2  yield  19.3e9  poisson  0.30 

mix  3 

ende 

* 

tracts 
pressure 
pfmix  -30.0e9 
pfvoid  -30.0e9 
pfracl  -35.0e9 
pfrac2  -20.0e9 
endf 


71 


Intentionally  left  blank. 


72 


APPENDIX  J: 

L/D  1/8  CTH  INPUT  DECK 


73 


Intentionally  left  blank. 


74 


*eor*  cgenin 
* 

cth  l/d=.125  (1/8  wa  impacting  rfia  at  2km/s 

control 

ep 

mmp 

*  usessd 
endcontroi 
convct 

convection  1 
interface  high 

endc 

* 

mesh 

block  1  geom  2dc  type  e 
xO  0. 

xl  n300  dxf  0.03175  rat  1.0 
endx 

yO  -.47625 

yl  n  320  dyf=.03175  rat  1.0 
endy 

xact  0.0  2.0 

yact  -.5  .5 

endbiock 

endmesh 

* 

insertion  of  material 
block  1 

package  'target' 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yl  0.0  y2  1000.0 

endinsert 

endpackage 

package  'rod' 

material  2 

numsub  25 

velocities  yvel  2.0e5 

insert  box 

xl  0.0  x2 1.27 

yl  -.3175  y2  0.0 

endinsert 

endpackage 

endbiock 

endinsert  ion 
* 

edit 
block  1 
expanded 


endblock 

endedit 


eo5  num  2 
matl  -1 
mat2  -2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1 .69  0. 
aneosi  3.1416  0.  -1.  -7.36e10  0.  0.  0. 
aneosi  0.  0. 0.  0.  0.  0.  0. 
aneosi  26  1 .0 

aneosw  -2  tungsten  rhug=-1  thug=-1 

aneosw  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneosw  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneosw  0.  0.  0.  0.  0.  0.  0. 

aneosw  74.  .9 

aneosw  28  .07 

aneosw  26  .03 

endeos 

* 

tracer 
block  1 

add  0.0-0.01 

add  0.0-0.3174 

add  0.0-0.15875 

endtracer 

* 

endinput 

* 

*eor*  cthin 

cth  l/d=.125  (1/8)  wa  impacting  rha  at  2km/s 
♦ 

control 

mmp 

tstop  10.0e-05 
cpshiti  999. 
rdumpf=14400 

*  pvoid  10.0 

*  usessd 

endcontro! 

* 

edit 

shortt 

time  0.  dtfrequency  1 .0e-06 
ends 
longt 

time  0.  dtfrequency  1 .0e-05 
endl 
plott 

time  0.  dtfrequency  1 .0e-05 
endp 
histt 
time  0. 


dtfrequency  0.05e-06 


time  2.0e-Q6  dtfrequency  0.1e-06 

time  10.0e-06  dtfrequency  1.0e-06 

htracerl 

htracer2 

htracer3 

endh 

ende 

* 

boundary 

bhydto 

block  1 

bxbot  0 

bxtop  1 

bybot 1 

bytop  1 

endb 

endh 

endb 

* 

epdata 

matep  1  yield  7.0e9  poisson  0.33 

matep  2  yield  19.3e9  poisson  0.30 

mix  3 

ende 

* 

fracts 

pressure 

pfmix  -30.0e9 

pfvoid  -30.0e9 

pfracl  -35.0e9 

pfrac2  -20.0e9 

endf 

* 
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Intentionally  left  blank. 


78 


APPENDIX  K: 


UD  1/16  CTH  iNPUTDECK 


79 


Intentionally  left  blank. 


% 


80 


*eor*  cgenin 
* 

cth  l/d=.0625  (1/16)  wa  impacting  rha  at  2km/s 
* 

control 

ep 

*  mmt 

*  usessd 

endcontrol 

* 

mesh 

block  1  geom  2dc  type  e 
xO  0. 

xl  n  300  dxf  0.015875  rat  1.0 

endx 

yO  -.3175 

yl  n  320  dy<=.015875  rat  1.0 
endy 

xact  0.0  2.0 

yact  -.5  .5 

endblock 

endmesh 

* 

insertion  of  material 
block  1 

package  ’target’ 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yl  0.0  y2  1000.0 

endinsert 

endpackage 

package  ’rod’ 

material  2 

numsub  25 

velocities  yvel  2.0e5 

insert  box 

xl  0.0  x2 1.27 

yl  -.15875  y2  0.0 

endinsert 

endpackage 

endblock 

endinsert  ion 
* 

edit 

block  1 

expanded 

endblock 

endedit 

* 

eos  num  2 
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matl  -1 
mat2  -2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1 .69  0. 
aneosi  3.1416  0.  -1.  -7.36e10  0.  0.  0. 
aneosi  0.  0.  0.  0.  0.  0. 0. 
aneosi  26  1.0 

aneosw  2  tungsten  rhug=-1  thug=-1 

aneosw  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneosw  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneosw  0.  0.  0.  0.  0.  0.  0. 

aneosw  74.  .9 

aneosw  28  .07 

aneosw  26  .03 

endeos 

* 

endinput 

* 

*eor*  cthin 

cth  i/d=.0625  (1/16)  wa  impacting  rha  at  2km/s 
* 

control 
tstop  8.0e-05 
cpshift  999. 
rdumpf=14400 

*  pvoid  10.0 

*  usessd 

endcontrol 

* 

edit 

shortt 

time  0.  dtfrequency  1 .0e-06 
ends 
longt 

timeO.  dtfrequency  1.0e-05 
endl 
plott 

time  0.  dtfrequency  1 ,0e-05 
endp 
histt 

time  0.  dtfrequency  1 .0e-06 

endh 

ende 

* 

boundary 
bhydro 
block  1 
bxbot  0 
bxtop  1 
bybot  1 
bytop  1 
endb 
endh 
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endb 


epdata 

matep  1  yield  7.0e9  poisson  0.33 

matep  2  yield  19.3e9  poisson  0.30 

mix  3 

ende 

* 

tracts 
pressure 
pfmix  -1.0e6 
ptvoid  -1.0e6 
pfracl  -20.0e9 
pfrac2  -30.0e9 
endf 


83 


Intentionally  left  blank. 


84 


APFtNDIX  L: 

L/D  1/32  CTH  INPUT  DECK 
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# 

*eor*  cgenin 
* 

cth  l/d=.03125  (1/32)  wa  impacting  rha  at  2krrt/s 
* 

control 

ep 

mmp 

*  usessd 
endcontro! 

mesh 

block  1  georri  2dc  typ  ■  e 
xO  0. 

xl  n  300  dxf  0.0079375  rat  1.0 
x2  n  75  ±\  0.0079375  rat  1.05 

end* 

yO  -.15875 

yi  n  320  dyf=. 0079375  rat  1.0 
y2  n  75  dyf=.0079375  rat  1.05 
andy 

xact  0.0  2.0 
yact  -.5  .5 
endblock 

endmesh 

* 

insertion  of  material 
block  1 

package  ’target’ 

material  1 

numsub  25 

insert  box 

xl  0.0  x2  1000.0 

yi  0.0  y2  1000.0 

endinsert 

endpackage 

package  ’rod' 

material  2 

numsub  25 

velocities  yvel  2.0e5 

insert  box 

xl  0.0  x2  1.27 

yi  -.079375  y2  0.0 

endinsert 

endpackage 

endblock 

endinseriion 

* 

edit 
block  1 
expanded 
endblock 
ended  it 
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eos  num  2 
mati  -1 
mat2  -2 

aneos  -1  iron  rhug=-1  thug=-1 
aneosi  1  -1  7.86  0.  0.  4.5e5  1 .69  0. 
aneosi  3.1416  0.  -1.  -7.36e10  0.  0.  0. 
aneosi  0.  0. 0.  0.  0.  0.  0. 
aneosi  26  1.0 

aneosw  -2  tungsten  rhug=-1  thug=-1 

aneosw  3  -1  17.3  0.  0.  4.0e5  1.43  0. 

aneosw  1.295  0.  -1.  -4.52e10  0.  0.  0. 

aneosw  0.  0.  0.  0.  0.  0.  0. 

aneosw  74.  .9 

aneosw  28  .07 

aneosw  26  .03 

endeos 

* 

tracer 

block  1 

add  0.0  0.0 

add  0.0  -0.079375 

add  0.0  -0.0396875 

endtracer 

* 

endinput 

* 

*eor*  cthin 

cth  l/d=.03125  (1/32)  wa  impacting  rha  at  2km/s 

* 

control 

mmp 

tstop  50.0e-06 
cpshift  999. 
rdumpf=14400 

*  pvoid  10.0 

*  usessd 
endcontroi 
convct 

convection  1 
interface  high 

endc 

* 

‘restart 

*  cycle=6644 

‘endrestart 

* 

edit 

shortt 

time  0.  dtfrequency  1 .0e-06 
ends 
longt 

time  0.  dtfrequency  1 ,0e-05 
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endl 

plott 

time  0.  dtfrequency  1 .0e-05 
endp 
histt 

time  0.  dtfrequency  0.05e-06 

time  2.0e-06  dtfrequency  0.1e-08 

time  10.0e-06  dtfrequency  1 .0e-06 

htracerl 

htracer2 

htracer3 

endh 

ende 

* 

boundary 
bhydno 
block  1 
bxbot  0 
bxtop  1 
bybot  1 
bytop  1 
endb 
endh 
endb 

epdata 

matep  1  yield  7.0e9  poisson  0.33 

matep  2  yield  19.3e9  poisson  0.30 

mix  3 

ende 

* 

tracts 
pressure 
pfmix  -30.0e9 
pfvoid  -30.0e9 
pfracl  -35.0e9 
pfrac2  -20.0e9 
endf 
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LIST  OF  SYMBOLS 


c  -  Material  Acoustic  Velocity 

cp  -  Penetrator  Material  Acoustic  Velocity 

ct  -  Target  Material  Acoustic  Velocity 

D  -  Penetrator  Diameter 

E  -  Penetrator  Kinetic  Energy 

H  -  Penetration  Channel  Diameter 

L  -  Penetrator  Length 

P  -  Penetration  Channel  Depth 

r  -  Target  radial  dimension 

t  -  Time 

v  -  Particle  Velocity 

Vs  -  Impact  Velocity 

p  -  Material  Density 

pp  -  Penetrator  Material  Density 

pt  -  Target  Material  Density 
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ATTN:  PLOTARD  Jean-Paul 
Boite  Postale  No.  7 
77161  Country 
TRANCE 

1  PRB  S.A. 

ATTN:  M.  Vansnick 
Avenue  de  Tervueren  168,  Bte.  7 
Brussels,  B-1150 
BELGIUM 


Organization 

AB  Bofors/Ammunition  Division 
ATTN:  Jan  Hasslid 

Anders  Nordell 
BOX  900 
S-691  80  Bofors 
SWEDEN 

Messerschmitt-Bdlkow-Blohm 
Dynamics  Division 
ATTN:  Manfred  Held 
P.O.  Box  1340 
D-8898  Schrobenhausen 
GERMANY 

TNO  Prins  Maurits  Laboratory 
ATTN:  H.  J.  Passman 
P.O.  Box  45 
2280  AA  Rijswijk 
Lange  Kleiweg  137 
THE  NETHERLANDS 

ISL 

ATTN:  CHANTERET  Pierre-Yve 
F68301  Saint-Louis 
C6dex,  12,  rue  de  I’lndustrie 
B.P.  301 
FRANCE 

Lothar  Meyer 
IFAM 

Lesumer  Heerstrasse  36 
Bremen  77 
GERMANY 
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